Supersolid in a one-dimensional optical lattice in the presence of a harmonic trap 
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We study a system of ultra-cold atoms possessing long range interaction (e.g. dipole-dipole in- 
teraction) in a one dimensional optical lattice in the presence of a confining harmonic trap. We 
have shown that for large enough on-site and nearest neighbor interaction a supersohd phase can 
be stabilized, consistent with the previous Quantum Monte Carlo and DMRG results for the ho- 
mogeneous system. Due to the external harmonic trap potential the supersolid phase coexists with 
other phases. We emphasize on the experimental signatures of the various ground state phases in 
the presence of a trap. 
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I. INTRODUCTION 



The realization of the supersolid form of matter, where 
the superfluid and the crystalHne order co-exist [l|, 0], 
in the ultra-cold bosonic atoms in optical lattices is at 
the forefront of research. After the claim of observing 
the supersolid phase in solid "^He by Kim ei a/ the 
progress in the research of this exotic phase of matter 
has advanced substantially. The successful observation 
of the superfluid (SF) to Mott insulator (MI) transition 
in ultra-cold bosonic atoms in 3-D 4] and subsequently 
m 2D ^\ and ID [6*] has shaped the study of ultra-cold 
systems as an ideal tool to understand condensed mat- 
ter phenomena. In order to achieve the supersolid form 
of matter that is characterized by the co-existence of the 
superfluid and crystalline order, it is essential for the sys- 
tem to have long range interactions. The remarkable ex- 
perimental realization of BEG in Cr atoms Q that have 
fairly large dipole moment, has increased the expecta- 
tions to observe the supersolid phase in optical lattice 
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experiments. 

In recent years there have been several theoretical evi- 
dences for the supersolid phase in various lattice geome- 
tries (sl-fT^. However, the experimental search of the su- 
persolid in the ultra-cold atomic systems in optical lat- 
tices still remains a challenge. The real experimental 
situation is different from the usual homogeneous system 
considered in theoretical calculations. In experiments, 
the translational symmetry of the lattice is broken due 
to the presence of an external harmonic trap potential 
(magnetic or optical) and various quantum phases co- 
exist [T5l - [2^ . Hence, it is essential to understand the 
signatures of the supersolid phase in the presence of such 
a trap. 

In this paper we have considered a system of ultra-cold 
bosonic atoms possessing long range interactions in a one 
dimensional optical lattice with a harmonic confinement. 
The Hamiltonian for this kind of system is represented 
by the extended Bose-Hubbard model, 
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<i,j> i 

+ V n^n-i +VtY^ rfn,. (1) 



<i,J> 



Here, t is the hopping amplitude between the nearest 
neighbor sites o|(ai) is the bosonic creation (an- 

nihilation) operator obeying the Bosonic commutation 
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relation [ui, Oj] = Sij and rii = ajui is the number oper- 
ator. U and V are the on-site and the nearest neighbor 
interactions, respectively. Vt is the magnitude of the ex- 
ternal trap potential and is the distance from the trap 
center. We re-scale in units of the hopping amplitude, i, 
setting t = 1, making the Hamiltonian and other quanti- 
ties dimcnsionless. 

The homogeneous version of this model (i.e. without 
the external trap), has been studied earlier using several 
techniques in one dimension m-Hii. The predic- 

tion of an accurate phase diagram using Quantum Monte 
Carlo method [s!] and DMRG [9| has revealed the phys- 
ical conditions required to stabilize a supersolid phase. 
It has been shown that the supersolid phase is obtained 
when: 

1. The total density of the system is incommensurate 
to the lattice. 

2. The on-site (U) and the nearest neighbor interac- 
tions {V) are fairly large compared to the hopping 
amplitude (t). 

3. The condition U < 2V is satisfied. 

A homogeneous system exhibits a uniform phase deter- 
mined by the global chemical potential for a given set of 
interaction parameters. The phase diagram of the model 
in Eq. [T]in the homogeneous limit i.e., Vt = 0, exhibiting 
different possible phases including the supersolid phase 
is shown in Fig. [1] Q ■ In the presence of an external trap, 
the role of a local chemical potential becomes important 
as demonstrated in our earlier work on the Bosc-Hubbard 
model [1^ . Since the local chemical potential varies from 
the center of the trap to the edges, the system exhibits 
different phases simultaneously. An earlier DMRG study 
of model given in Eq.[T] could not confirm the presence of 
the supersolid phase in the system ,32t | . A recent study of 
this model in two dimension using mean field theory pre- 
dicts that the noise correlation could be a valid signature 
to separate the supersolid phase from the other ground 
state phases [ll]. In this paper we re- visit the extended 
Bose-Hubbard model with the external harmonic trap 
potential and search for experimental signatures of the 
different ground state phases, in particular, the super- 
sohd phase. 

The remaining part of the paper is organized as follows. 
In Sec. ini we discuss the method of our calculation using 
the finite size density matrix renormalization group (FS- 
DMRG) technique. The results along with discussions 
are presented in Sec. IIIII with experimental signatures 
for the different ground state phases in Sec. IIVI and we 
present our conclusions in Sec. [V] 

II. METHOD OF CALCULATION 

To obtain the ground state of model ([1]) for the system 
of N bosons on a lattice of length L, we use the FS- 
DMRG method with open boundary conditions [s^ [sj] . 
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FIG. 1: (Color on line) The phase diagram of the model in 
Eq.[T]in the homogeneous limit i.e., Vt = and for U = 10.0 
in the \i — V plane. 
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FIG. 2: (Color on line) Homogeneous phase diagram for the 
model in Eq. [T] showing canonical trajectories. For a given 
value of V, the presence of an external trap allows all phases 
that fall on the line that starts from the canonical trajectory 
to the V-axis for example: lines AB and CD in the figure. 

This method has been widely used to study the Bose- 
Hubbard model m, H, m Hi, m. we have considered 
six bosonic states per site and the weights of the states 
neglected in the density matrix formed for the left or the 
right blocks are less than 10~^ [1^. In order to improve 
the convergence of the results, the finite-size sweeping 
procedure as given in [2^ [s^l has been used for every 
length. Using the ground state wave function \iIiln) and 
energy El{N), we calculate the following physical quan- 
tities and use them to identify the different phases. 
The on-site local number density (rii), defined as, 

(rii) = (■i/'L^v|f^i|^L^f), (2) 
gives the local density distribution. The fluctuation in 
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the local number density, k^, which is finite for the SF 
phase, is calculated using the relation 

«i = - {"n-if (3) 

and finally the existence of the CDW order is confirmed 
by calculating the structure factor: 

S{k)^j-^Y.^''^'-^Hn,n,). (4) 

In our calculations, we have considered a system of length 
L — 140 and vary TV from 30 to 140. In our previ- 
ous work on the homogeneous extended Bose-Hubbard 
Model [q^ , we had considered a fixed value of the on-site 
interaction U = 10 and vary the nearest neighbor inter- 
action strengths V from to 10. We choose the same 
range of parameters here as well, since the homogeneous 
phase diagram for this range, as shown in the Fig. [U ex- 
hibits most of the interesting phases for this model. The 
strength of the external confining trap potential is fixed 
at Vt = 0.008. 



III. RESULTS AND DISCUSSION 

We begin with the summary of the phase diagram for 
the homogeneous extended Bose-Hubbard model, which 
has been studied recently [1, for a wide range of den- 
sities and interaction parameters namely, the on-site in- 
teraction U and the nearest neighbor interaction V . The 
phase diagram for a typical value of the on-site interac- 
tion, say U = 10 is shown in Fig. [1] Q . The phase dia- 
gram consists of gapped as well as gapless phases. The 
gapless phases include the superfluid phase, the super- 
solid phase where superfluidity and charge density wave 
order co-exist, and the solitonic phases. The gapped 
phases are (i) the Mott insulator phase with p = \ 
for V < Vc 5.4, (ii) charge density wave phase 
CDW-II, (where every other site is doubly occupied, i.e., 
|2 2 • • •)) with average density p = 1 for F > Vb ~ 5.4 
and (iii) the CDW-I phase (alternative sites are occupied, 
i.e., boson density varies as |1 1 • • •)) with average 
density p = 1/2 for V > Vc ^ 3.0. The gap vanishes 
when doping above or below these gapped phases. For 
example doping below half-filling (p = 1/2) gives rise to 
solitons that break the CDW-I order. This phase ex- 
tends over a small range of densities below the CDW-I 
and eventually goes over to the superfluid phase when 
the density is further decreased. However, the behavior 
of the system when doping above half- filling is different. 
For small V we get similar solitonic phases, however, 
for larger V a supersolid phase stablizes. The super- 
solid phase forms again while doping above and below 
the CDW-II phase. In fact there exists a range densities 
0.5 < p < 1 and p > 1 for V > U/2 where the supersolid 
phase is the stable ground state of model ([T]) as shown in 
the Fig. m 



Let us introduce a harmonic trap potential. Earlier 
studies of the one-dimensional Bose-Hubbard model in 
the presence of a trap has demonstrated the co-existence 
of the superfiuid and the Mott insulator phases [l^, [13] . 
The Mott insulator is characterized by the formation of a 
plateau in the local number density (n^) as a function of 
the distance from the center of the trap and is incom- 
pressible, while the superfiuid phase is characterized by 
large local number density fluctuations and is compress- 
ible. The nearest neighbor interaction brings about the 
charge density wave order in the system due to the in- 
terplay between the U and V terms in the Hamiltonian. 
Figures [3] and m show the density profile, i.e., the varia- 
tion of the local density (n;) as a function of the distance 
from the trap center r^. We obtain the density profile 
for two sets of parameters: (i) for the number of bosons 
fixed at N ^ 80, but different values of V (Fig. ^ and 
(ii) fixed nearest neighbor interaction, V = 8, but differ- 
ent values of N (Fig. 2]). The following three features are 
clearly seen: (i) the local density (rii) is maximum at the 
center of the trap, (ii) the density falls-off with increase 
in r.i and (iii) the density profile exhibits plateaus and 
oscillations. 

In order to understand these features and identify var- 
ious phases from the density profile, we define the local 
chemical potential at the site i at a distance r; from the 
center of the trap as. 

Mi = Mo - Vrr-f. (5) 

Here = El{N -I- 1) — El{N) is the chemical potential 
of the system. For the homogeneous system, Hi = 
for any i. However, for a finite trap the local chemical 
potential pi equals moi which is its maximum value, at 
the center of the trap and decreases radially outward as 
in Eq. O It is instructive to plot the density profile as a 
function of pi instead of r.i as in Fig. [5l It may be noted 
from Fig. [5] that the density of bosons at any site i is 
controlled by the value of the local chemical potential pi. 
So a decrease in pi results in a decrease in (n;), with the 
maximum at the center of the trap as observed in Figs. [3] 
andm 

The density of bosons plays a very crucial role in the 
determination of the ground state of the model in Eq. [T] 
The gapped phases are possible only when the density is 
commensurate. The homogeneous system with a given 
value of U and V and a uniform local chemical potential 
Pa represents one point in the phase diagram. However, 
for the system with a trap potential, the density varies 
across the lattice due to the variation of the local chem- 
ical potential and therefore different phases co-exist. In 
order to understand this feature of co-existence of the 
different ground state phases and the role played by the 
local chemical potential, we study first the path in the 
phase diagram that is traced by pq as we change the 
interaction parameter V keeping the number of bosons 
A^, the trap potential Vt and the on-site interaction U 
fixed. This path is referred to as the Canonical Trajec- 
tory since N is held fixed. Fig. [5] shows several 



4 




r, 
1 



FIG. 3: The local density {jii) as a function of the distance 
from the center of the trap for iV = 80, (7 = 10, Vr = 0.008, 
but for different values of V . 
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FIG. 4: The local density {rn) as a function of the distance 
from the center of the trap for fixed 1/ = 8, f/ = 10, 
Ft = 0.008, but for different values of A^. 



canonical trajectories (for different values of N) in the 
homogeneous phase diagram. In may be noted that /iq 
is the local chemical potential at the center of the trap 
and the position of these canonical trajectories trace the 
phase present at the trap center as V is varied for fixed 
N . For example when N = 30, the canonical trajectory 
and hence the phase at the center of the trap goes from 
the superfluid to CDW-I as we increase V . The posi- 
tion of the canonical trajectory in the phase diagram can 
be shifted by changing the number of bosons TV. When 
the number of bosons is increased, say to iV = 40, /io in- 
creases and the position of the canonical trajectory in the 
phase diagram is shifted upward. As a result the center 
of the trap, say for y = 0, which was in the SF phase for 
= 30, is now in the Mott insulator phase. Following 
the canonical trajectory for N = 40, the trap center goes 
from MI to SF and then to a supersolid phase for increas- 



ing V . Thus the position of the canonical trajectory for 
a given N and V in the phase diagram represents the 
phase at the center of the trap. 

Moving away from the center of the trap, the local 
chemical potential decreases as in Eq. ([5]) and the varia- 
tion of /ii is represented in the phase diagram by a line 
drawn vertically downwards from the canonical trajec- 
tory to the horizontal axis. The local chemical potential 
values across the lattice fall on this line, which passes 
through different ground state phases. Therefore, the 
local chemical potential (and thus local density) at dif- 
ferent sites favor the co-existence of different phases in 
the presence of a trap. It is useful to re-plot the density 
profile given in the Fig. [3] as a function of /i^ using Eq. [5] 
instead of as in Fig. [S] We also calculate and plot, in 
the same figure, the average local number density define 
as 

fij = ((2nj -I- rii+i -f n^_i))/4. (6) 

For iV = 80 and V — 2.0, falls in the superfluid 
phase above the p — 1 Mott lobe (point A as indicated 
on the canonical trajectory corresponding to TV = 80 
in Fig [2) . This means that the center of the trap has 
(rii) > 1. Moving away from the trap center, /Zi de- 
creases along the line AB and there are regions where /i^ 
falls inside the MI lobe. From Fig.[5l we see that for these 
values of fj,i, {rii) = 1. Similarly as we move towards the 
edge, the values of decreases further such that the sys- 
tem is once again in a superfluid phase on the lower side 
of the Mott lobe. So the system for TV = 80, I^ = 2 has a 
superfluid core flanked by a MI phase and finally ending 
with a superfluid edge. The density profile (top panel of 
Fig [3] and Fig[5|) correlates with this result. In addition, 
there are oscillations in {rii) in the superfluid shoulders 
near fii — 1/2. The reasons for these oscillation are the 
following. For V = 2, the system is close to CDW-I lobe 
(see Fig[T]). In the thermodynamic limit, the CDW-I or- 
der can stabilize only for V > Vc ~ 3.0. However, the 
finite size of the system allows a CDW-I phase to exist 
for lower values of V, here V = 2, although it vanishes 
in the thermodynamic limit. Finite size effects are char- 
acterized by oscillations in the local density (rii). These 
oscillations stabilize at higher values of V into the CDW- 
I phase. For example kr V = 7 and N = 80, fiQ (point C 
in Fig.HD falls inside the CDW-II lobe yielding a CDW-II 
phase at the center. As we move towards the edges, the 
CDW-II phase is flanked by a supersolid phase, CDW-I 
and finally a superfluid shoulder as can also be infered 
from the density proflles as shown in Fig. [3] and Fig[S] In 
fact, these conclusions can be further fortified by com- 
paring the variation of the average as a function of 
fii with the density of the corresponding homogeneous 
system as in Fig. [51 The agreement is striking, leading 
to the conclusion that for a given set of parameters, the 
phase of a system with an external trap is represented 
by a line starting from the canonical trajectory to the 
horizontal axis while the phase of the homogeneous sys- 
tem is represented by a point in the phase diagram. This 
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immediately shows that while the homogeneous system 
can have a unique phase, the phases tend to co-exist for 
an inhomogeneous system. 

In addition to scanning along the phase diagram at 
fixed values of N and varying V , it is also equally pos- 
sible to fix the nearest neighbor interaction V and move 
along the phase diagram by varying N and therefore the 
chemical potential /iq. The canonical trajectory in the 
phase diagram moves upwards (downwards) by increas- 
ing (decreasing) the total number of bosons and as a re- 
sult the local chemical potential at the center of the trap 
fiQ changes, giving rise to different phases at the center. 
This is demonstrated in Fig. 2] for fixed y = 8 for differ- 
ent values of N. For N — 30, position of the /io is inside 
the CDW-I lobe (see Fig. and as discussed above, the 
corresponding system has a CDW-I core flanked by a su- 
perfluid edge as seen in the density profile (top panel 
of Fig. SI). An interesting situation occurs for N — 40, 
where the trap center is expected to be in the elusive su- 
persolid phase, as seen in Fig. [5] and is characterized by 
density fluctuations between 1.0 < {rii) < 1.5, that is, 
the system has a CDW order at incommensurate densi- 
ties [9] . As a result, the system now will have a supersolid 
core, followed by a CDW-I and a superfiuid phase mov- 
ing outward from the trap center (top panel of Fig. |4|) . 
Further increase in N leads to the inclusion of a CDW- 
II phase in the system in addition to the supersolid, the 
CDW-I and the superfiuid phases. 
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FIG. 6: (Color on-line) Average local number density for 
system with a trap and (rii) for a homogeneous system as 
function of the local chemical potential fj,i for V — 2 and 7. 




FIG. 7: (Color on-line) Number density per site and its fluc- 
tuations that serve as a measure of compressibility and hence 
can be used as a tool to pick out the compressible and the in- 
compressible phases that coexist in the presence of a harmonic 
trap. 



FIG. 5: (Color on-line) Local number density {rii) and average 
define as rii = {{2ni + n^+i -I- ni_i)/4) as a function of local 
chemical potential for different values of V but fixed A'' = 80. 

The next issue we address here is a scheme to pick out 
the various phases using local properties of the system. 
We will follow the discussions in '2^ and use local com- 
pressibility or equivalently the fluctuations in the number 
density per lattice site, Ki, given in Eq. [31 as a tool to 
distinguish between the gapped and the gapless phases. 
It is known that the number fluctuation is large in the 
superfluid phase while it is a minimum for the MI and 
the CDW phases. Fig. [7] shows the variation of Ki across 



the lattice. For small values of V, Ki varies at the cen- 
ter and at the edges of the trap indicating that these 
regions are in the superfluid phase, while the plateaus 
represent the Mott insulator phase. Further, we note 
that these plateaus (minima) occur exactly over the val- 
ues of ri where the average local density Ui exhibits a 
plateau at integer densities. Therefore, one can pick out 
the incompressible phases using the density proflle and 
its local fluctuation Ki and identify them using the phase 
diagram and the canonical trajectories. As an example. 
Fig. |8] shows the different phases for iV = 80 but varying 
V. In the next section, we will discuss the experimental 
signatures for the various phases that have been isolated 



SF ; SS 
CDW-I 



CDW - II 



SS ; SF 
CDW-I 



-0.5 



-60 



-40 



-20 



20 



40 



60 



0.2 



0.5 1 



1.5 2 2.5 0.5 
P 



1.5 2 2.5 0.5 
P 



1.5 2 2.5 0.5 
P 



1.5 2 2.5 
P 



FIG. 8: (Color on-line) Picking out the different phases from 
the density profile for iV = 80 and V = 8.0. 

in the presence of a harmonic trap using global properties 
of the system. 



IV. EXPERIMENTAL SIGNATURES 

The presence of a harmonic trap in the optical lattice 
leads to the co-existence of the superfluid, the Mott insu- 
lator, the charge density wave and the supersolid phases 
as seen in the previous sections. As a result, extract- 
ing the signature of a particular phase in the presence 
of other phases becomes a theoretically important exer- 
cise in order to make connections with experiments. In 
the following we analyze possible global signatures of the 
various ground state phases that can be experimentally 
confirmed. 

It is now possible in experiments to record the spatial 
distribution of the lattice with different filling factors [35- 
[38| . Similar experiments in one-dimensional optical lat- 
tices can yield density profiles using which the ground 
state phases can be mapped. Another way to obtain di- 
rect information about the Mott plateaus (shells in 3D) 
is through the atomic clock shift experiment [s^. By 
using density dependent transition frequency shifts, sites 
with different occupation can be spectroscopically distin- 
guished, thus giving us information about the number of 
sites corresponding to a given density p of bosons, de- 
fined as N{p). As a first step, we look for the signatures 
of the solid phases (MI and CDW) in an atomic clock 
shift experiment. 

In Figs. l9l and [TOl we plot N{p) as a function of p for 
different values of V fixing = 80 and different N val- 
ues with fixed V = 8 respectively. The density profiles 
corresponding to these parameter values are given in the 
Figs. [3] and |4] respectively. The presence of the incom- 
pressible phases, that is, the MI, the CDW-I and II in the 
system can be inferred from the formation of a peak in 



FIG. 9: N{p) versus p for A' = 80 and different values of the 
nearest neighbor interaction V. The presence of incompress- 
ible phases can be distinguished by the formation of a peak 
at commensurate densities. 
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FIG. 10: N{p) versus p for V — 8.0 as number of Bosons 
is varied. Incompressible phases can be picked out by the 
formation of peaks at commensurate densities. 



N{p) at commensurate densities. For example, existence 
of a Mott plateau in the density profile for V ranging be- 
tween and 5 ( see Fig[3l) correlates with a peak in N{p) 
at p — 1. Similarly peaks in N(p) at p = 2 correlate with 
the formation of CDW-II phases in the density profile. 
Similar conclusions can be drawn from Fig 1101 Compar- 
ing with the density profile in Fig|4l we can conclude that 
the formation of peaks in N(p) at integer densities can 
be correlated with the existence of the solid phases, i.e., 
MI or CDW. 

In order to distinguish between the two solid phases, 
i.e, the CDW and MI phase, we calculate the structure 
factor, as defined in Eq. |4l Fig. [TT] shows the structure 
factor in momentum space as V is varied for N = 40, 
while Fig. [T^] has fixed V = 8.0 for different N values. 
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FIG. 11: Structure factor as a function of q for = 40 and 
different values of V . 



that can be used to check for the presence of a SS phase 
in the trap is the momentum distribution n{q) (ioj . 
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FIG. 12: Structure factor as a function of q for V = 8.0 and 
different values of A*'. 



FIG. 13: (Color on-line) Momentum distribution as a function 
of q for A'^ = 40 and different values of V. 
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From the phase diagram (see the canonical trajectory in 
Fig. ^ ioi N — 40 the phases at low V values are the su- 
perfluid and the Mott insulator. However, for higher val- 
ues of V, a CDW-I phase is possible. The CDW oscilla- 
tions in the density profile translates to the formation of a 
peak at q = TT in the structure factor. As V increases, this 
peak at g = TT grows in magnitude reaching its maximum 
value when the trap center exhibits the CDW crystalline 
structure. However this crystalline structure is possible 
for a CDW or a SS phase. For example, for TV = 40 
and V = 8.0, the center of the trap is in the supersolid 
phase that has the CDW-like crystalline structure and 
is compressible like a superfluid. Hence the next step 
is to distinguish between the CDW ordered phases that 
could be either compressible (SS phase) or incompress- 
ible (CDW phase itself). While this can be established 
locally with the behavior of compressibility as a function 
of the distance from the trap center, a global signature 



FIG. 14: (Color on-line) Momentum distribution as a function 
of q for V = 8.0 and different values of A^. 

In experiments, the bosons in the optical lattice are 
allowed to expand and the interference pattern in the 
density is recorded. The density distribution is mirrored 
in the momentum distribution defined as, 

1 ^ 

^ lYI exp{iq{k - 0) (7) 

kj=l 

which then provides global information about the various 
phases present in the system. Figures [T51 and [HI show the 
momentum distribution, respectively, for various values 
of V with fixed = 40 and for various values of N with 
fixed V = 8.0. We see that in addition to the peaks at 
q = and q = ±27r, there are peaks around q = n. In or- 
der to understand the reason for this peak at g = tt, let us 
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look at the momentuni distribution for the homogeneous 
system as in Fig. [151 We choose four densities to demon- 
strate the features of the peak in n{q) at q = tt. From the 
phase diagram we see that for V = 8.0 the homogeneous 
system with p = 0.42 is in the superfluid phase, p = 1/2 
and 1 are, respectively, in CDW-I and CDW-II phases 
and for p — 0.67, the system is in the supersolid phase. 
From Fig. [15] we note that the presence of a supersolid 
order in the system is accompanied by a peak in the mo- 
mentum distribution at q = tt, which is absent in the SF, 
CDW-I and CDW-II phases. The structure function for 
the same set of densities as given in Fig. [THj show peak 
at q = IT when the system is in CDW-I, CDW-II and SS 
phases. This confirm that the peak in n{q) at q = vr is a 
clear signature of the supersolid phase. 

Therefore for the trapped systems, when the phases 
co-exist, we note that a peak in the momentum distribu- 
tion function at g = tt signals the presence of a supersolid 
phase somewhere in the trap, and the peak height being 
maximum when the supersolid occupies the center of the 
trap. We summarize below the signatures of the differ- 
ent ground state phases for the inhomogeneous extended 
Bose-Hubbard model: 

• MI - Phase: 

— Peaks in N{p) at integer densities. 

— No peaks in the momentum distribution at 
n{q — tt). 

— No peaks in the structure function at S{q = 
n). 

• CDW Phase: 

— Peaks in N{p) at integer densities. 

— No peaks in the momentum distribution at 
n{q ~ vr). 

— Peaks in the structure function at S{q = tt). 

• supersolid phase: 

— No Peaks in N{p) at integer densities. 

— Peaks in the momentuni distribution at n{q = 
n) 

— Peaks in the structure function at S{q = tt) 



the system. We show the signature of different phases by 
calculating different observable quantities such as the on- 
site number density, the number fluctuation, the struc- 
ture factor and the momentum distribution. We also 
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FIG. 15: (Color on line) Momentum distribution for homoge- 
neous case. Note that when the system is in the Supersolid 
phase, a peak a.t q — tt develops. 




0.5 
0.4 
^ 0.3 
^ 0.2 
0.1 




1 , 1 

1 1 p = 0.67 


I ' 1 

1 


I i 

- 1 , ! , J 


1 i- 



FIG. 16: (Color on line) Structure function for homogeneous 
case. Note that when the system is in the Supersolid phase, 
a peak at q — n develops. 



V. CONCLUSION 

In conclusion, we have studied a system of dipolar 
ultra-cold bosonic atoms in the frame work of the ex- 
tended Bose-Hubbard model in the presence of external 
harmonic trap. Using finite size density matrix renor- 
malization group (FS-DMRG) method we have demon- 
strated the simultaneous existence of different phases in 



document global signatures for the ground phases that 
can be observed experimentally. 
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